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The reduced cross section for exclusive (e, e'p) reactions has been studied in DWIA for the ex- 
ample of the nucleus 16 using a spectral function containing effects of correlations. The spectral 
£^ , function is evaluated directly for the finite nucleus starting from a realistic nucleon-nucleon inter- 

action within the framework of the Green's function approach. The emphasis is focused on the 
correlations induced by excitation modes at low energies described within a model-space of shell- 
model configurations including states up to the sdg shell. Cross sections for the p-wave quasi-hole 
■ transitions at low missing energies are presented and compared with the most recent experimental 

data. In the case of the so-called perpendicular kinematics the reduced cross section derived in 
DWIA shows an enhancement at high missing momenta as compared to the PWIA result. Further- 
more the cross sections for the s- and d-wave quasi-hole transitions are presented and compared to 
available data at low missing momenta. Also in these cases, which cannot be described in a model 
without correlations, a good agreement with the experiment is obtained. 
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I. INTRODUCTION 



o 
o 

The quasi-elastic (e, e'p) reaction constitutes a very well suited tool to study the limitations of the simple shell- 
model or independent-particle model (IPM) of the nucleus. Indeed, it is now generally accepted that atomic nuclei are 
*G ' many-body systems in which correlations beyond the mean-field or Hartree-Fock picture play a significant role. It has 
been argued that the strong short-range and tensor components of realistic nucleon-nucleon (NN) interactions should 
induce short-range correlations into the nuclear wave function. These correlations should give rise to an enhancement 
in the momentum distribution of quasi-hole states at high momenta as compared to the momentum distribution 
£h ■ derived from a Hartree-Fock or mean field description of the nucleus. Therefore high-resolution (e, e'p) experiments 
have been performed to determine the spectral function of nucleons at high momenta leading to the ground-state or 
• i-h , states with low excitation energies in the daughter nucleus |I]J^] . 

Microscopic calculations, which account for the effects of these short-range correlations, indeed predict components 
in the momentum distribution at momenta around 2-3 fin -1 , which are larger by orders of magnitude as compared to 
the predictions of Hartree-Fock or IPM calculations. These high-momentum components, however, mainly originate 
from the spectral function at large missing energies, leading to configurations of the (A-l) particle system at energies 
above the threshold of particle emission for another nucleon. The momentum distribution for nucleons at low missing 
energies, which can be explored in exclusive (e, e'p) reactions, seems not to be very sensitive to short-range correlations 
and may rather well be approximated by those derived from the IPM |3]-||]. Similar results are also obtained in the 
study of the spectral function for nuclear matter [pHl2"[ . 

This analysis was essentially confirmed by experiment. The experimental data for the reduced cross section in 
the light nucleus 16 at low missing energies could be described reasonable well within the prediction of an IPM 
[jij. On the other hand, the momentum distribution for the heavy nucleus 208 Pb showed a small enhancement of 
high momemtum components at low missing energies as compared to a Hartree-Fock or IPM prediction This 
enhancement is not as large as the one predicted from short-range correlations for the total momentum distribution 
including the contributions from high missing energies. It has been argued that this deviation from the IPM might 
be due to long-range correlations, which corresponds to low-energy excitations of the many-body system P,p~3| . 
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The fact that this enhancement of the momentum distribution is observed for a heavy nucleus but not for the light 
nucleus 16 supports the idea that the enhancement at small energies may originate from long-range correlations. 
The effects of short-range correlations should be rather independent of the nuclear system under consideration as 
these correlations are not very sensitive to the global structure of the whole nuclear system. Consequently the effects 
of short-range correlations should be rather similar for the nuclei 16 and 208 Pb. Contrary long-range correlations 
could be sensitive to the whole nuclear system and exhibit different results for different nuclei. They are related 
to the excitation modes at low energy and therefore results derived from nuclear matter, which shows a continuous 
single-particle spectrum, can be quite different from those in finite nuclei, for which the low-energy excitations are 
rather sensitive to the shell-structure. Guided by these considerations we investigated long-range correlation effects in 
the spectral function and other related observables directly for the finite nucleus under consideration jl4j , as a study 
of nuclear matter may not be very reliable. 

If one wants to study the effects of correlations by a comparison with experimental data, one has to employ the 
correlated wave function in a complete description of the (e, e'p) reaction. Going beyond the Plane Wave Impuls 
Approximation the cross section cannot be factorized any longer into the elementary electron-proton cross section 
and the spectral function, exhibiting the effects of correlations. The reduced cross section calculated in the Distorted 
Wave Impuls Apoproximation (DWIA) may contain effects of the final state interactions (FSI) and other effects which 
cannot be separated in this approach from the effects of correlations. Therefore a careful analysis has to be made of 
the different components entering in the computation of the reduced cross section. 

Also we want to emphasize that one should investigate the reduced cross section leading to final states in the 
daughter nucleus, which are absent in the IPM. In the case of 16 processes of this kind would be the knock-out of 
a nucleon leading to the 5/2 and 1/2 states of positive parity with low excitation energy in 15 N. In the IPM the d 5 / 2 
and lsj/2 states are not occupied and therefore the corresponing momentum distributions are not dominated by the 
quasi-hole part of the simple shell-model. 

To account for long-range correlations a harmonic oscillator basis is used in order to determine the energies and 
the mixing of shell-model configurations. We consider that the low-energy excitation modes are adequately described 
within a model-space which includes all orbitals up to the sdg shell. A finite basis of oscillator states, however, is not 
at all appropriate to describe high-momentum components in the nuclear wave function, since these high-momentum 
components will be dominated by the tail of the oscillator basis states. Therefore, as we will explain below, we have 
used a mixed representation of basis states, which considers a shell-model basis to describe the excitation modes, 
but a basis of plane-wave states to determine the spectral function. The Green's function formalism for the nuclear 
many-body theory (for a reference see e.g. the recent review articles [ p^|p^ |) will be applied to evaluate the spectral 
function. 

In Sec. [n] we will present the main points of the method we use to compute the correlated spectral function for 
16 0. This part is followed in Sec. Ill by a brief summary of the formalism used to describe the (e,e'p) reaction, 
where Coulomb distortion effects of the electron wave function and final state interaction of the outgoing proton are 
considered. The results for the reduced cross section for the different approximations and partial wave transitions are 
presented in Sec. |j"v| . At the end we give a short review of our main conclusions. 



II. THE SPECTRAL FUNCTION 



The calculation of the cross section for exclusive (e, e'p) reactions requires the knowledge of the hole spectral 
function. In the case of finite system it is convenient to introduce a partial wave decomposition which yields the 
spectral function for a nucleon in the single-particle basis with orbital angular momentum I, total angular momentum 
j, isospin t and momentum k 

S ljT (k,k';E) = Yl (*ol4'yl*£ _1 >(*£~V«r|*£> S(E - (E* - E^)) , (1) 

n 

where akij T ( a \nj T ) denotes the corresponding annihilation (creation) operator. The state Y&q) refers to the ground 
state of the target nucleus, while l 1 !^ -1 ) is used to identify the nth excited eigenstate of the hamiltonian with one 
particle removed from the target nucleus. Hence the hole spectral function in its diagonal form 5(k, E) gives the 
probability of removing a particle with momentum k from the target system of A particles leaving the resulting (A-1) 
system with an energy E A ~ l = Eq — E, where Eq is the ground state energy of the target. The spectral function for 
the various partial waves, Sij T (k, k' , E) can be obtained from the imaginary part of the corresponding single-particle 
Green's function gij(k\ 1 k 2 , E). Note that here and in the following we have dropped the isospin quantum number r, 
as we ignore the Coulomb interaction between the protons and study a symmetric nucleus with N = Z. 
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To determine the correlated single-particle Green's function one has to solve a Dyson equation, which using the 
partial wave representation appropriate for finite systems can be written as 



gij(ki, ki; w) = g\j '(h, k 2 ;u) + / dk a J dk^g^ '(ki, fc 3 ; u)AT,i j (k 3 , fc 4 ; w)gy(fc 4 , k 2 ; u) , (2) 

where g^ HF ' refers to the Hartree-Fock propagator and AEy represents contributions to the irreducible self-energy, 
which go beyond the Hartree-Fock approximation for the nucleon self-energy used to derive g^ HF ^ . The hole spectral 
function Sij can then be calculated easily from the imaginary part of the single-particle Green's function by 

Sij(k,k',u>) = -Imag[#y a;)] . (3) 

Although the evaluation of the Hartree-Fock approximation to the Green's function g[f F ^ (k\ , fe; oS) , the definition of 
AS;, and the technique used to solve the Dyson Eq.(|^) have been already discussed in detail in previous publication 
0,0, we include a brief summary of the relevant aspects of this method. 

A. Model space and effective hamiltonian 

Long range correlations are taken into account by means of the Green's function approach within a finite model 
space. This model space shall be defined in terms of shell-model configurations including oscillator single-particle 
states up to the sdg shell. The oscillator parameter, b — 1.76 fm, has been chosen appropriate for the nucleus 16 0. 
This model space does not allow the description of short-range correlations. Nevertheless, we also have to take into 
account the effects of short-range correlations by introducing an effective interaction, i.e. a G- matrix appropriate for 
the model space. This truncation of the Hilbert space into a model space, the degrees of which are treated explicitly, 
and the space outside this model space, which is taken into account by means of effective operators, is often referred 
to as a double partitioned Hilbert space and has been used before for finite nuclei Jl8],[l9| and nuclear matter po| , pT| . 

The G-matrix is determined as the solution of the Bethe-Goldstone equation 

Q = V + V _ n Qmo * G, (4) 

w ^mod 1 ^mod 

where T is identified with the kinetic energy operator, while V stands for the bare two-body interaction. For the 
latter we have chosen the Reid soft-core potential |^2| . In this equation the Pauli operator Q moc j is defined in terms 
of our harmonic oscillator single-particle states. Thus applying Q moc j to two-particle states \a(3 > one obtains 

!0 if a or (3 below Fermi level 

if a and (3 in model space (5) 
| af3 > elsewhere 

The model space used in the Eq. ^ and ||) includes all states up to the sdg-shell. Note that with this definition 
of Q moc i we ensure that no doublecounting of correlations occurs between the short-range correlations taken into 
account in the ^-matrix and the long-range correlations evaluated by means of the Green's function method within 
the model-space. 

In the solution of the Bethe-Goldstone Eq.(Q) we have chosen a constant value of lu — —30 MeV for the starting 
energy. This value is a reasonable mean value for the sum of two single-particle energies for hole states in 16 O. Clearly, 
this choice is an approximation introduced to simplify the calculations, but one has to note that our results do not 
depend significantly on the actual value of u. The use of a constant starting energy also implies that we do not try to 
account for a depletion of the occupation probability due to scattering into states outside the model space as it has 
been done e.g. in p^p^ |. 

The matrix elements for Q are computed in a plane-wave basis for a specific finite nucleus and a correspondig 
mode l-sp ace by expanding the techniques for solving the Bethe-Goldstone equation for finite systems as described in 
Ref. || . 
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B. Nucleon Self-Energy and Green's Function 



The calculation of the self-energy is performed in terms of two-particle states, characterized by single-particle 
momenta in the laboratory frame. Such a antisymmetrized 2-particle state would be described by quantum numbers 
such as 

\kxhhk 2 hhJT > (6) 

where ki , U and ji refer to momentum and angular momenta of particle i whereas J and T define the total angular 
momentum and isospin of the two-particle state. The transformation from the relative and cm. coordinates, in which 
the matrix elements of Q are defined to the states displayed in Eq.@ can be made by use of the well known vector 
bracket transformation coefficients J2(||27j . 

Performing an integration over one of the ki, one obtains a 2-particle state in a mixed representation of one particle 
in a bound harmonic-oscillator state while the other is in a plane-wave state 

oo 

\n1hj1k2hj2 JT > = J dki k\ R 7ll h{ bk i) \k1hj1k2hj2 JT > . (7) 



Here R ni h stands for the radial oscillator function. An oscillator length b = 1.76 fm ( Uw = 13.3 MeV) has been 
selected, which is an appropriate value to describe the bound single-particle states in 16 O. Now with the help of eqs.(^ 
- 0) we can write down our Hartree-Fock approximation for the self-energy in momentum representation 

^ 31 ' n 2 l 232 JT 

The summation over the oscillator quantum numbers is restricted to the states occupied in the IPM of 16 0. This 
Hartree-Fock part of the self-energy is real and does not depend on energy. One obtains the HF single-particle wave 
functions by expanding them 

\a HF ljm >= \Kiljm >< K i \a HF >y (9) 

i 

in a complete and orthonormal set of regular basis functions within a spherical box of radius Rbox which is large 
compared to the radius of the nucleus 

* Ujm (r) = (r\Kiljm) = N aJl {K ir )y ljm {'dip) (10) 

where Nu is an appropriate normalization constant, yij m denotes the spherical harmonics including the spin degrees 
of freedom while ji stands for the spherical Bessel functions with discrete momenta Ki determined from the boundary 
condition 

3i(K t R box ) = 0. (11) 
The expansion coefficients of Eq. (g) are obtained by diagonalizing the HF Hamiltonian 

E <**l + ( K n\a H % = e*J (Kila^. . (12) 

n— 1 

Here and in the following the set of basis states in the box has been truncated by assuming an appropriate N. From 
the HF wave functions and energies one can construct the HF approximation to the single-particle Green's function 
in the box, which comes out as 

ffi ( »,y iB) . <*lgf^<p*>it . (13) 

Note that by choosing especially this basis we are able to separate contributions from different momenta to the HF 
single-particle state, which is essential in order to compute at the end of our formalism the single-particle Green's 
function in momentum space. 
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FIG. 1. Contributions to the nucleon self-energy of second (a) and third order in the interaction (b) and (c) 

The next step in the evaluation of the irreducible self-energy is to take into account terms of second order in the 
effective interaction which correspond to intermediate 2-particle 1-hole (2plh) states as displayed in Figjlja) 



E (2plh),, ,/ s _ 1 \ - < kh\Q\p x p 2 >< PiP2\G\k'h > 

13 l ' ,UJ> ~9. ^ ^ u-e(p 1 ,p 2 ,h) + iri 



h<F p ly p 2 >F 

and intermediate 2-hole 1-particle (2hlp) states 



p>F h!,h2<F 



< k'p\Q\h 1 h 2 >< h 1 h 2 \Q\k''p > 
lu - e(hi,h 2 ,p) - ir\ 



Here we have introduced the abbreviation 



e(a,/3, 7 )=er +e 



HF , HF HF 



(14) 



(15) 



(16) 



where the are the HF single-particle energies. Note that the summations in Eq.( |lj, |T^) on particle labels like 

Pi , P2 and p are restricted to those single- particle states within the model space, which are above the Fermi level (F) , 
whereas the labels hi, h 2 and h refer to hole states. 

After the definition of the self-energy we can now proceed and calculate the corresponding single-particle Green's 
function gij by solving a Dyson equation ( see Eq.(|^)) with the g HF taken from Eq.(|l3|) and including the correlation 
effects contained in ASjj 



(17) 



C. Solution of the Dyson Equation 



The technique we use to solve the Dyson equation in order to extract the basic ingredients of the single-particle 
Green's function is very similar to the one developed in |l^,|l^,Q. So we will restrict ourselves in giving a short 
review of the basic steps towards the determination of the single-particle Green's function for a finite system within 
a model space of discrete single-particle states. 

In order to obtain the information necessary for the Lehmann representation of the Green's function, we rewrite 
the Dyson equation as an eigenvalue problem [fL9[ 



an 
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An 

\Ain 
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V Y Q J 



(18) 
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where for simplicity we have dropped the corresponding conserved quantum numbers for parity and angular momen- 
tum (Ij). The matrix to be diagonalized contains the HF hamiltonian defined in ( |l2| ) and the coupling to the P 
different 2p\h configurations and Q 2h\p states which can be constructed in our model space with quantum numbers 
for parity and angular momentum j, which are compatible to the single-particle quantum numbers Ij under consider- 
ation. As long as we are still ignoring any residual interaction between the various 2p\h and 2hlp configurations the 
corresponding parts of the matrix in ( |l8| ) are diagonal with elements defined by e; (Ej) for 2p\h (2hlp) 

ei = e(p!,p2,h) Ej = e(7ii,/i 2 ,p) , (19) 

where we have used again the abbreviation (|l6| ) . The matrix elements connecting the HF part to the additional states 
refer to 

a m i =< k m h\g\pip 2 > 

A m j =< k m p\Q\hih 2 > (20) 

Solving the eigenvalue problem (Eq.([l8])) one gets as a result the single-particle Green's function in the Lehmann 
representation in the discrete basis of the box defined in Eq. (|Io|) . The eigenvalues ui n define the position of the poles 
of the Green's function, which refer to the various states of the system with A±l nucleons and the corresponding 
spectroscopic amplitudes are given by 

<*oKK +1 >=Xo, ki for u n >E F 

<H\*l i \<- 1 >=Xlk i for oj n <E F (21) 

which depend on whether u n is an energy above or below the Fermi energy Ep. Note that the coefficients Xq k . in 
the above equations stand for the momentum representation of the quasi-hole (quasi-particle) wave functions. That 
means one can set X£ k . = & n ij{ki). With the help of this nomenclature we can finally write down the spectral 
function in momentum space for a given energy w„ and a given partial wave 

Sy (fc, k', lu u ) = X^ k , Xl k = cj>* nlj {k') 0„y (k) (22) 

in a separable representation (compare Eq. ([!])), which is important in order to be able to use the so computed spectral 



function in the description of (e, e'p) reactions. The corresponding spectroscopic factor S n (u> n ) (compare Eq.(29)) for 
the removal of a particle of a given shell {nlj} is determined by the norm of the quasi-hole wave function which reads 
according to the chosen box basis as 

SnM = f drr 2 |$„y(r)| 2 = £ \X^ k f , (23) 
o 

where $ n y(r) stands for the Fourier-Bessel transform of <fr n ij(k). 

In a straightforward way one can improve the approximation discussed so far and incorporate the effects of residual 
interactions between the 2p\h configurations as illustrated in the diagrams displaying the self energy in Fig. [j}}) and 
c). The same holds for the 2h\p configurations. One simply has to modify the corresponding parts of the matrix in 
Eq.(|l8|) and replace 

'e x ... \ /Ei ... \ 

; ' • . J U 2p ih , and : J H 2 hi P , (24) 

... e P ) V ... Eq) 

where H.2 P ih and HihXp contain the residual interactions in the 2p\h and 2h\p subspaces. The solution of the 
eigenvalue problem also leads to a normalization condition, which ensures that 

E \ x ^f + E raj 2 = E i < > i 2 + E i < ^ki*? > i 2 = i . (25) 

n m n m 

where the sum on n accounts for all solutions with uj n larger than the Fermi energy and the sum on m for all solutions 
with u n below the Fermi energy. Again this implies that one ignores all effects of correlations, which are due to 
configurations outside the model space, like e.g. an effective energy- dependent hamiltonian |23]p4|| . In that case one 
has to renormalize the condition of Eq. ( p5| ) as well. 

Note that for the solution of the eigenvalue problem one can apply the so-called "BAsis GEnerated by Lanczos" 
(BAGEL) scheme |l9| , |i"7| , p8| in order to get a very efficient representation of the single-particle Green's function in 
terms of a few "characteristic" poles in the Lehmann representation. 
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III. REDUCED CROSS-SECTIONS 



Reduced cross-sections for (e, e'p) reactions on 16 have recently been measured at NIKHEFK |2^] in the region of 
missing momenta up to around 300 MeV and at MAMI in Mainz jij including missing momenta up to 700 MeV. The 
simplest approximation to analyze the (e, e'p) process is the Plane Wave Impulse Approximation (PWIA), where one 
makes the assumption that the proton is ejected from the nucleus without any further interaction with the residual 
nucleus. In nonrelativistic PWIA the differential cross section factorizes into two terms, the elementary electron- 
proton cross section, accounting for the interaction between the incident electron and the bound proton, and the 
spectral function that accounts for the probability to find a proton with given energy and momentum in the nucleus. 
Although the factorization is destroyed when one takes into account the distortion of the electron and/or outgoing 
proton waves, or a relativistic approach for the bound state, it is useful and common practice to analyze the results in 
terms of a reduced cross section defined in such a way that it would coincide with the spectral function if factorization 
were fulfilled. For selected values of the missing energy E m (i.e. for selected quasi-hole excitations) the reduced cross 
section is given by 

= jL dEm W^' 1 dE Jj^ e ' < 26 > 

with p m the experimental missing momentum (= —Pa—x), E p , |p p |,f2p (e' e ,Q' e ) the outgoing proton (electron) kine- 
matical variables, and experimentally, the integral is performed over the interval AE m that contains the peak of the 
transition under study. The term a ep represents the elementary electron-proton cross section. The data of p(p m ) 
are obtained dividing the experimental cross section by o^, as given by Eq.(17) of Ref. J|(|. We therefore use the 
same expression for a ep in our theoretical calculations. In PWIA, p(p m ) represents the momentum distribution of 
the selected quasi-hole state a. 

In this section we briefly summarize the formalism used to describe the (e, e'p) reaction. More details can be found 
in Refs. |3l|-|33]|. We base our calculations on the impulse approximation (virtual photon absorbed by the detected 
nucleoli), which is known p3] to be a reliable approximation at quasi-elastic kinematics. The basic equations that 
determine the reduced cross section are given explicitly in Refs. in terms of the electron and nuclear currents. 

The calculations have been performed with the code developed by one of us J35| . 

We work in the laboratory frame in which the target nucleus is at rest and use the notation and conventions of 
Ref. [Q. We denote by kg = (e e ,k e ) the four-momentum of the incoming electron and by = (e' e ,k' e ) the four- 
momentum of the outgoing one. The four-momentum of the exchanged photon is q^ L — kg—k"^ = (w,q). P% = (Ma,0) 
and P^-i — [E'A-ii'P A-i) denote the four-momenta of the target and residual nucleus, while p£ = (E p ,p p ) is the 
four-momentum of the ejected proton. 

Using plane waves for the electrons and considering knock-out from a given {nlj} quasi-hole state, we write the 
amplitude for the (e, e'p) process in DWIA as |3^-|34|]: 

Wij = -^=u(k' e ,a' e )^u(k e ,a e )tllr N{ujj(i ) , (27) 

where u(k, er) represent four-component relativistic free electron spinors and J^(w,q) is the nuclear current 

J£(w,q)=j;^ 5(E F -Er-u>) fdre*' <^ F A \J^ A > , (28) 

IF J 

where the matrix element of the nuclear charge-current density operator is taken between the initial > and the 
final > nuclear states. The initial state will be given by the many-particle wave function of A bound particles in 
the ground state for the target nucleus \^>q >. For the final state, the experimental conditions dictate a state that 
behaves asymptotically as a knocked out nucleon with momentum p„ and a residual nucleus in a well-defined state 
l*^" 1 ^) > with energy E and quantum numbers n. It is possible [pp7|] to describe the matrix element in Eq. ( J2q ) 
in terms of the simple one-body current operator sandwiched between the hole spectral functions 

[s n }h En (p) =< *H-He)Mp)\H > (29) 

and 

X p ; ) En =<*t 1 (EMp)\<> (30) 
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describing the overlap between the residual state \'^^~ 1 (E) > and the hole produced in >, and |\&^ > respectively, 
by removing a particle with momentum p. The way to compute these overlaps was described in the previous section. 
Hence, we take 4>En normalized to 1, and S n (E) is the spectroscopic factor associated with the removal process. 

One may perform a similar calculation for \ p e > the nna l (continuum) wave funciton of the ejected proton. Usually, 
however, this particle wave function is derived from a phenomenological local optical potential. Also in this work 
we compute the the wave function for the outgoing proton Xp e as a scattering solution obtained from an optical 

potential in a relativistic framework fitted to elastic proton scattering data on 16 O pafa . No Perey factors are included, 
as the dependence on the energy of these potentials is very soft at the energies of interest. 

To avoid expansions in p/M and the reductions of the current operator typical of the nonrelativistic approaches, 
the fully relativistic kinematics and structure of the operator is used. In order to do that, we build a 4-component 
spinor out of the nonrelativistic bispinor (f>E n - Also, we use relativistic optical potentials |38| ], which according to 
previous results [ j32|]39| , seem to be more adequate for (e, e'p) calculations. However, our calculation is essentially 
nonrelativistic and, in consistency with the derivation of the spectral functions, no contributions from the negative 
energy sector are allowed in neither the bound or scattered proton wave function. The results would be equivalent 
to the ones obtained with the nonrelativistic formalism and equivalent optical potentials, only that all orders of p/M 
are included in our calculation. We obtain similar results with a nonrelativistic prescription for the scattered states 
(nonrelativistic optical potentials and Perey factors), only the reduced cross-sections are slightly increased. 

For the nucleon current operator we take the free nucleon expression 

J» = F,r + i-^^Q* > (31) 

where F\ and F^ are the nucleon form factors related in the usual way J3(| to the electric and magnetic Sachs form 
factors of the dipole form. As discussed in Refs. [p^2pc| l, DWIA results depend on the choice of the nucleon current 
operator. Here we have chosen the operator that is closer to the one used in the nonrelativistic calculations of the 
dweepy code, and most often employed. 

The Coulomb distortion of the electron wave function is considered with an effective momentum approach. The 
value of the shift in the energy of the electron is obtained by comparing with the exact calculation (DWBA) |3j| and 
it is set in this case to about 4 MeV. For a light nuclei as 16 O, this procedure is already sufficient in order to include 
this effect accurately in the cross-sections ED. 



IV. RESULTS AND DISCUSSION 



In the first part of this discussion we would like to emphasize the effects of the final state interaction for the 
outgoing proton. For this purpose Fig. |^ displays the reduced cross section, calculated in perpendicular kinematics for 
the 16 0(e, e'p) 15 N reaction leading to the ground state (1/2 - ) of 15 N, as a function of missing momentum. Results 
are shown for three different momentum transfers q of the virtual photon. As a reference all three parts of this figure 
include the results of the PWIA. The PWIA yields a reduced cross section which is identical to the spectral function 
and independent on the momentum of the absorbed photon. Two different models have been employed to describe 
the distortion of the wave function of the outgoing proton. The first model is based on a microscopic Dirac-Brueckner- 
Hartree-Fock calculation for the nucleus 16 O |)42|| . This DWIA-RH (Relativistic Hartree) approach only yields a real 
part for the nucleon nucleus potential. The second approach (DWIA-ROP) also accounts for the absorption of the 
outgoing proton by means of inelastic scattering. Also the DWIA-ROP approach assumes a relativistic form of the 
otical potential. The scalar and vector terms for the real and imaginary part are obtained in a phenomenological way 
by a global fit of elastic proton- nucleus scattering data |38| . 

With the inclusion of the final state interaction, the cross section cannot be factorized any more into the free electron- 
nucleon cross section times a spectral function. Therefore the reduced cross sections displayed in Fig. 2 depend on 
the momentum q of the absorbed photon. The difference between the two-models of the final state interaction can 
mainly be attributed to the effects of the imaginary component in DWIA-ROP. The absorptive effects contained in 
this imaginary part leads to a reduction of the cross section (as compared to the DWIA-RH results) which seems 
to be rather insensitive with respect to the missing momentum p m . Therefore in the following we will restrict the 
presentation of calculated cross sections to the DWIA-ROP approach. 

The real part of the optical potential tends to smoothen the reduced cross section by shifting strength from small 
momenta to high missing momenta. Independent on the detailed form of the optical potential one also observes an 
increase of the reduced cross section with increasing photon momentum q. 
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FIG. 2. Reduced cross section for the ie O(e, e'p) 15 N reaction leading to the ground state of 15 N in perpendicular kinematics. 
The results for PWIA are compared with those from DWIA calculated with a relativistic Hartree potential (RH) and a 
relativistic optical potential (ROP). 

After exploring these features of the final state interaction, we now want to compare our results for the reduced 
cross sections with experimental data. In particular we would like to study how the effects of long-range correlations 
may effect the shape of the spectral function. The long-range correlations are described in terms of the 2plh and 
2hlp configurations within a model space including single-particle states up to the sdg shell and taking into account 
effects of the residual interaction between these configurations. For a more detailed discussion on the importance 
of contributions to the self-energy represented by diagrams of third and higher order in the residual interaction as 
displayed by Figs. [l]b and [l]c and their effects on the momentum distribution and other observables we refer to Ref. 

To observe the influence of long-range correlations on the momentum distribution, we present in Fig. rathe reduced 
cross section for the 16 0(e, e'p) 15 N reaction leading to the ground state (l/2~) and first excited state (3/2 - ) of 15 N. 
These have been calculated both with the fully correlated spectral function (Full) and also with our mean-field (HF) 
description (see Eq.(p"3|)). The results are compared with the experimental data taken from the experiment at MAMI 
(Mainz) J5J. The underlying formalism to compute the reduced cross section is performed as described previously 
where the phenomenological relativistic optical potential was used. The relevant piece of the nuclear structure 
calculation for the proton knock-out reaction is the quasi-hole part of the fully correlated spectral function for the 
Pi 12 and P3/2 partial wave presented in the previous chapter (see Eq.(p2|) and Eq.(|24|)). As we are mainly interested 
in the comparison of the data at high missing momenta with the predictions at low momenta, the corresponding 
spectroscopic factors S n (E) 0.60 (Opi/2) and 0.45 (0p 3 / 2 ) were determined by a fit of the calculated cross section to 
the experimental results at small missing momenta deduced from the NIKHEF data [^9|. These adjusted spectroscopic 
factors are significantly smaller than the theoretical values of 0.83 and 0.85 for the p 1 / 2 and p 3 / 2 states, respectively. 
For this comparison, however, one must keep in mind that these calculated spectroscopic factors only account for the 
effects of long-range correlations. Short-range correlations should lead to additional reduction by another factor of 
0.8 and subnucleonic degrees of freedom may be responsible for the remaining discrepancy. 
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FIG. 3. Reduced cross section for the 0(e,e'p) N reaction leading to the ground state (l/2~) and first excited state 
(3/2~) of 15 N in the kinematical conditions considered in the experiment at MAMI (Mainz) [l|. Results for the mean-field 
description (HF) and the fully correlated spectral function (Full) are presented. The spectroscopic factors determined by a fit 
to low-p data from the NIKHEF experiment [^9) are S , o Pl/2 =0.60 (0.83) and 5o P3/2 =0.45 (0.85), where the values obtained by 
the theoretical approach are enclosed in parenthesis. 

As one can clearly see in Fig. ^, the inclusion of long-range correlations does not lead to a significant enhancement 
of the reduced cross section at high missing momenta as compared to the predictions of the mean field or Hartree-Fock 
approach. This conclusion is not altered even in a complete description of the (e,e'p) experiment including FSI and 
other effects. This statement confirms the result already obtained by an investigation of long-range correlations (i.e. 
correlations induced by excitation modes at low energies), on the energy and momentum distribution for the nucleus 
16 [jl4j. It should be emphasized, however, that a rather reasonable description of the reduced cross section is 
obtained over a large range of missing momenta with the adjustment of only one parameter, the spectroscopic factor, 
which has been adjusted to describe the NIKHEF data at small missing momenta. 

Now one could argue that in the (3/2~) and (l/2~) case the spectral function has a dominant quasi-hole part 
originating from the corresponding hole (i.e. bound) state. Therefore any effects of correlations would be overwhelmed 
by this dominating quasi-hole component, which is present already in the Hartree-Fock approximation. This argument 
is supported by observing that in Fig. || the reduced cross section computed with the fully correlated spectral function 
does not differ substantially from the mean field description. In this sense it would be more interesting to investigate 
transitions to final states which do not have such a dominant quasi-hole part and would be impossible to describe 
within the mean field approximation. Best examples at low missing energies in 16 is the case of the c? 5 / 2 and S1/2 
hole states. As we are able to compute the fully correlated spectral functions also for these partial waves it would be 
worth comparing these results with corresponding experimental data. Unfortunately, until now we are not aware of 
data from the MAMI (Mainz) experiment Q for high missing momenta of these partial waves. So we had to take the 
experimental data from the NIKHEF experiment in parallel kinematics [ p9| which cover only the momentum region 
up to 280 MeV/c. Moreover, they could not resolve the contributions from the (5/2 + ) and (l/2 + ) states. In Fig. [| we 
display the results obtained for the reduced cross section taking the fully correlated spectral functions and assuming 
an incoherent sum of the l/2 + and 5/2 + contributions for the 16 0(e, e'p) 15 N reaction leading to the first excited 
state with positive parity (l/2 + — 5/2 + -dublett) of 15 N. The experimental data points are taken from the NIKHEF 



experiment |29| mentioned above. 

Contrary to the previous case, the curves in Fig. || are scaled with the spectroscopic factors given by the theoretical 
approach, 0.055 (si/2) and 0.035 (^5/2) respectively. The reasonable agreement with experiment shown in Fig. [|, means 
that the underlying microscopic calculation of the spectral function could be used to explore the high momentum 
region. Note that one could get better agreement with the experimental results by using the spectroscopic factors of 
0.0357 (si/2 and 0.1140 (d 5 / 2 which have been obtained by fitting the experimental data |29[] . In view of our aim to 
test our theoretical approach this has not been done here. 
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FIG. 4. Reduced cross section for the ls O(e, e'p) 15 N reaction leading to the first excited state with positive parity 
(l/2 + — 5/2 + -dublett) of 15 N in parallel kinematics. The experimental data are taken from NIKHEF experiment 0|. The 
corresponding spectroscopic factors Sa ± , 2 =0.055 and Sd 5/2 =0.035 are taken from the theoretical approach. 



The reduced cross section for the (e, e'p) process has been studied in DWIA for the example of the nucleus 16 using 
a spectral function containing long-range correlations. The fully correlated spectral function in the various partial 
waves is derived from the single-particle Green's function, which is obtained from the solution of a Dyson equation 
in basis of plane-wave states. The self-energy for the nucleons entering this Dyson equation contains the Brueckner- 
Hartree-Fock term plus the coupling to 2plh and 2hlp configurations inside a model space with single-particle states 
up to the sdg shell. The long-range correlations are described in terms of these 2plh and 2hlp configurations including 
the residual interaction among them. 

Effects of final state interactions for the outgoing protons are included by means of relativistic optical model 
interactions. The effects of the final state interaction are non-negligible and spoil the factorization of the total cross 
section into the spectral function times the free electron-nucleon cross section. The effects of the final state interaction 
turn out to be rather insensitive to the specific choice of the optical model. 

Including these effects of the final state interaction the calculated cross section for (e,e'p) reactions leading to the 
Opi/2 and 0p 3 / 2 quasi-hole states in 16 O agree rather well with the experimental data fit] over a large range of missing 
momenta. It turn out, however, that these data are rather insensitive to correlation effects. Spectral functions derived 
from the Hartree-Fock approach yield an agreement of similar quality. 

In order to investigate the effects of correlations, one should study (e, e'p) reactions to final states, which are 
impossible within the mean field or Hartree-Fock approach. Examples for such states at low missing energies are the 
5/2+ and l/2 + states in 15 N. The cross section derived from the correlated single-particle Greens function is in good 
agreement with the experimental data. 
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